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Thermodynamic properties of quasiparticles in a graphene-based structures are investigated. Two 
graphene superconducting layers (one superconducting component is placed on the top layered- 
graphene structure and the other component in the bottom) separated by oxide dielectric layers 
and one normal graphene layer in the middle. The quasiparticle flow emerged due to external 
gate voltage, we considered it as a gas of electron-hole pairs whose components belong to different 
layers. This is a striking result in view of the complexity of these systems: we have established that 
specific heat exhibits universal (-T 3 ) behaviour at low T, independent from the gate voltage and 
the superconducting gap. The experimental observation of this theoretical prediction would be an 
important step towards our understanding of critical massless matter. 
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INTRODUCTION 



Graphene is a unique system in many ways [1H3J. It 
is truly two-dimensional systems, has unusual electronic 
excitations described in terms of Dirac fermions that 
move in a curved space, is an interesting mix of a semi- 
conductor (zero density of states) and a metal (gapless- 
ness), and has properties of soft matter. The electrons in 
graphene seem to be almost insensitive to disorder and 
electron-electron interactions and have very long mean 
free paths 0]. Hence, graphene's properties are different 
from what is found in usual metals and semiconductors 
0, ||| . Graphene has also a robust but flexible structure 
with unusual phonon modes that do not exist in ordi- 
nary three-dimensional solids. In some sense, graphene 
brings together issues in quantum gravity and particle 
physics, and also from soft and hard condensed matter. 
Interestingly enough, these properties can be easily mod- 
ified with the application of electric and magnetic fields, 
addition of layers, control of its geometry, and chemical 
doping @, H|. Moreover, graphene can be directly and 
relatively easily probed by various scanning probe tech- 
niques from mesoscopic down to atomic scales, because 
it is not buried inside a 3D structure @. This makes 
graphene one of the most versatile systems in condensed- 
matter research. Besides the unusual basic properties, 
graphene has the potential for a large number of ap- 
plications [l(| , from chemical sensors to transistors 
[illil. Graphene can be chemically and/or structurally 
modified in order to change its functionality and hence- 
forth its potential applications. Moreover, graphene can 
be easily obtained from graphite, a material that is abun- 
dant on the Earths surface. This particular characteristic 
makes graphene one of the most readily available mate- 
rials for basic research. Whereas many papers have been 
written on monolayer graphene in the past few years, only 



Vtg 



_3fc— -„.**■- - 



Va 



Vb 



* e-mail: d.bolmatov@qmul.ac.uk 



Vbg 



FIG. f: Graphene multilayer device tunneling structure. Two 
sheets of superconducting graphene are separated by thin di- 
electric oxide layers and one normal graphene monolayer in 
the middle. Separate gate electrodes make it possible to vary 
independently the carrier concentration in the normal and 
superconducting graphene layers. 



a small fraction actually deal with multilayers [lj, [l5| . 
The majority of the theoretical and experimental efforts 
have concentrated on the single layer, perhaps because of 
its simplicity and the natural attraction that a one atom 
thick material, which can be produced by simple methods 
in almost any laboratory, creates. Nevertheless, few-layer 
graphene is equally interesting and unusual with a tech- 
nological potential, perhaps larger than the single layer 
flU]. Indeed, the theoretical understanding and experi- 
mental exploration of multilayers is far behind the single 
layer. 

Graphene can be considered as a semiconductor with 
zero band gap. Electron energy spectrum of graphene 
contains two Dirac points that separate the electron and 
the hole subband. In a multilayer structure the Fermi 
levels of the hybrid system layers can be adjusted inde- 
pendently by the gate voltage. The electron- hole sym- 
metry near the Dirac points ensures perfect nesting be- 
tween the electron and the hole Fermi surfaces. A flow of 
electron-hole pairs in the graphene-layered structures is 



equivalent to two oppositely directed electrical currents 
in the layers. Therefore, the flow of such pairs is a kind of 
superconductivity [l7j . It is believed that electron- hole 
pairs may demonstrate superfluid behaviour [IH [l9j . In 
this letter we consider multilayer graphene structure and 
claim: flow of electron-hole pairs consisting of an elec- 
tron from the top (bottom) graphene layer and a hole 
from the bottom (top) graphene layer behave as super- 
conductive one. This phenomenon can be obtained in 
a graphene-layered structure with dielectric layers, pair- 
ing of electrons in one top (bottom) graphene layer with 
holes in the bottom (top) graphene layer occurs. 

In the weak-coupling limit, exciton condensation is a 
consequence of the Cooper instability [2(| of solids with 
occupied conduction-band states and empty valence- 
band states inside identical Fermi surfaces. Each layer 
has two Dirac-cone bands centered at inequivalent points 
in its Brillouin zone. The particle-hole symmetry of the 
Dirac equation ensures perfect nesting: the nesting con- 
dition requires only that the Fermi surfaces be identi- 
cal in area and shape and not that layers have aligned 
honeycomb lattices and hence aligned Brillouin zones, be- 
tween the electron Fermi spheres in one layer and its hole 
counterparts in the opposite layer, thereby driving the 
Cooper instability. Global wave vector mismatches can 
be removed by gauge transformations. When weak inter- 
valley electron-electron scattering processes are included 
only simultaneous momentum shifts of both valleys in a 
layer are allowed. 



II. GRAPHENE-BASED STRUCTURES 

Multilayered structures are the building blocks of many 
of the most advanced devices presently being developed 
and produced. They are essential elements of the highest- 
performance optical sources and detectors, and are being 
employed increasingly in high-speed and high-frequency 
digital and analog devices. The usefulness of such struc- 
tures is that they offer precise control over the states and 
motions of charge carriers. 

In Fig [1] we illustrate a device geometry in which flow 
of electron-hole pairs may be observed and in Figj2] we 
sketchily represent the dispersion relation of elementary 
excitations in graphene multilayer device tunneling struc- 
ture. We have two superconducting graphene monolayers 
which are separated by thin oxide barriers and one nor- 
mal graphene monolayer in the middle. Top and bottom 
gates (Vt g and Vb g ) are used to electrostatically. The top 
and bottom gates are separated manipulate the quasipar- 
ticle concentrations in the top and bottom layers. The 
top and bottom gates are separated from the graphene 
layers by gate oxides which are several nanometers thick. 
As the phenomenon of superfluidity of electron- hole pairs 
has not yet been experimentally observed in graphene 
multilayer structures. 

A dissipationless flow of electron-hole pairs in equilib- 
rium through the graphene mono-layer and thin oxide 



layers, depending on the phase difference between the 
two superconducting graphene layers. The single-particle 
Hamiltonian in graphene is the two-dimensional Dirac 
Hamiltonian 



H = 



H+ 
H_ 



where H± — —ihvF(a x d x ± cr y d y + U) acting on a four- 
dimensional spinor (^a+^b+^a-^b-)- The indices 
A, B label the two sublattices of the honeycomb lattice of 
carbon atoms, while the indices ± label the two valleys 
of the band structure. There is an additional spin degree 
of freedom, which plays no role here. The 2x2 Pauli 
matrices <7j act on the sublattice index. 

The time-reversal operator interchanges the valleys 



T = 



a z 
a z 



P = T- 1 



with p the operator of complex conjugation. We consider 
a sheet of graphene in the x — y plane. Electron and hole 
excitations are described by the Bogoliubov-de Gcnncs 
equation 



H - E F A 
A* E F - THT- 1 



= e 



with u and v the electron and hole wave functions, 
e > the excitation energy (relative to the Fermi en- 
ergy Ep), H the single-particle Hamiltonian, and T the 
time-reversal operator. The pair potential A(r) couples 
time-reversed electron and hole states. In the absence of 
a magnetic field, the Hamiltonian is time-reversal invari- 
ant, THT^ 1 = H and we yield two decoupled sets of 
four equations each, of the form 



H± — Ep A 
A* E F - H± 



= £ 



Separate gate electrodes make it possible to vary inde- 
pendently the carrier concentration of electron-hole pairs 
in the normal direction (z-direction) to graphene mono- 
layer, thin oxide layers and superconducting graphene 
layers. For — d < z < d, the pair potential vanishes 
identically, disregarding any intrinsic superconductivity 
of graphene, where d is a total width of graphene mono- 
layer and thin oxide layers. For z < —d and z > d the su- 
perconducting graphene layers will induce a nonzero pair 
potential A(z) via the proximity effect similarly to what 
happens in a planar junction between a two-dimensional 
electron gas and a superconductor [12] ■ The bulk value 
A e^ for z < -d and A e"^ for z > d (with <j> the 
superconducting phase) is reached at a distance from 
the interface which becomes negligibly small if the Fermi 
wavelength X F in superconducting layers is much smaller 
than the value X F in graphene mono- layer and thin ox- 
ide layers. We assume that the electrostatic potential U 
in graphene mono-layer, thin oxide and superconducting 
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FIG. 2: Scheme of dispersion relation of elementary exci- 
tations in graphene multilayer device tunneling structure. 
The spectrum of elementary excitations of quasiparticles in 
graphene-based structure was proposed in the He 4 -like man- 
ner: the linear portion near k — represents "phonons" (due 
to single normal graphene layer in the middle of structure) 
and the portion near k = ko corresponds to "rotons", which 
requires a minimal energy Ao for its creation. 



graphene layers may be adjusted independently by a gate 
voltage or by doping. Since the zero of potential is arbi- 
trary we may take: U(r) = for — d < z < d and U(r) — 
—Uq otherwise. For Uq large positive, and Ep > 0, the 



Fermi wave vector k. 



2tt/X f = (E F + U )/hv F 



graphene superconducting layers is large compared to the 
value Uf = 2h/\f = Ef/Hvf in graphene mono-layer 
and thin oxide layers (with vf the energy- independent 
velocity in graphene). 



III. SPECIFIC HEAT 

For the graphene mono-layer and oxide thin layers 
(where A = = U) we assume that Fermi level is tuned 
to the point of zero carrier concentration. In the super- 
conducting graphene layers there is a gap in the spectrum 
of magnitude |A| = A . 

The canonical partition function for this graphene- 
based structure can be interpreted as a collection of nor- 
mal mode oscillators, the oscillator labelled by k con- 
taining rik quanta of energy e(k). First we calculate the 
Helmholtz free energy for normal graphene mono-layer 
(phonon-like dependence) and thin oxide layers (its make 
dispersion relation between linear and quadratic depen- 
dence smooth, see Fig[2]) which can be written as 



F n (T, V) = -2Nk B TV 



dk, 
— I 
2tt 



-A; In 



1 + exp 



~k^r) 



Riemann zeta function. Contribution to the specific heat 
per unit volume is 



(1) 



where kd {—kd) is the largest (lowest) wavevector for 
which the linear dispersion is a reasonable approxima- 
tion. Therefore, the contribution to free energy from the 
normal region of graphene-based structure is 



F n (T, V) 



3((3)Nk B T 3 V 



2irh 2 



(2) 



where N is the number of the 2 component Dirac fla- 
vors, N = 4 in the single layer graphene and C( n ) is the 
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Q T 2 
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(3) 



A. Specific heat of superconducting layers: 
graphene-based layers 

The grand potential of the system (roton-like depen- 
dence) is given by 



q(V, T) ee -pF s (V, T) = — In [1 - exp (-/3e(k))] 

k 

~5>xp(-/3e(k)) ~iV (4) 

k 

where N is the "equilibrium" number of rotons in the 
superconducting graphene layers. The summation over 
k may be replaced by an integration, with the result 



F s = -Ak B TN 
= -8k B TV 



k 2 



dk 

2^ 



ke -f3{&+U+hv F {k-k )) ^ 



where k\ and k% are the lowest and highest values re- 
spectively of k for which the quadratic approximation to 
the dispersion curve is reasonable (the limits of integra- 
tion may be extended to ±oo). The free energy of the 
roton-like supcrfhud gas is given by 



F, 



4(k B T) 2 V (k B T 



ko 



-(2A cos 4>+U)/k B T 



(G) 



The roton-like contribution to the specific heat we esti- 



mate as: C- 



V 



V.l 



Cy 2 , where 



C 



4k B (k B T) 2 (2Ao< 



6 + 4- 



tt(Hvf) 2 
2 A cos <p + U 



i<j>+U)/k B T 



k B T 



2 A cos 4> + U 



_ Ak B k {k B T) _ (2Ao 

°V,2 — Z e 

irnvF 



2 + 2- 



2A cos 4> + U 
kWT 



k B T 

cos4>+U)/k B T x 
2 



/ 2A cos 4> + U 
{ k^f 



(7) 



(8) 



where fco = n/d, d — 0.4 nm is the thickness of graphene 
sheet, n is the number of oxide dielectric layers (the thick- 
ness of each of them is presumed equivalent to graphene 
one) between normal and superconducting graphene lay- 
ers. Rough low-temperature estimates often quote the 
result 



S cx CI 



CI 



(9) 



where S is thermopower [21| . In FigJ3]is depicted the spe- 
cific heat against temperature. It is remarkable that in 
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FIG. 3: Dependence of specific heat versus temperature ex- 
hibits Debye-like behaviour (T 3 -like) . This is a striking result 
in view of all four curves sit on top of each other. 



IV. CONCLUSION 

The practical significance of this investigation rests 
on the expectation that high-quality contacts between 
a superconducting graphene and normal graphene sheets 
can be realized. This expectation is supported by the 
experience with carbon nanotubes (rolled up sheets of 
graphene), which have been contacted succesfully by 
superconducting electrodes. Graphene-based structures 
provides a unique opportunity to explore the physics of 
the "relativistic Josephson effect", which had remained 
unexplored in earlier work on relativistic effects in high- 
temperature and heavy-fermion superconductors. 

We believe that this device structure can also lead to 
greater understanding in another condensed matter sys- 
tems such as high-temperature superconductors and ex- 
citon gases. The low cost experimental techniques allow 
for easy verification of the proposed hypothesis. 



graphene-based superconducting structure the phonon- 
assisted drag effect is suppressed. The carrier concentra- 
tion can be tuned varying different parameters: the thick- 
ness of oxide layers, phase difference in superconducting 
layers, pair and electrostatic potentials on one hand and 
on the another proposed structure exhibits Debye-like be- 
haviour. 
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